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ABSTRACT 


4 method is described for making estimates oi the total 
emissivity of hydrogen in the temperature and pressure ranges where 
hydrogen atoms predominate under equilibrium conditions. Fora 
characteristic geometrical depth of 50 cm, and temperatures of the 
order of 12, 500°K and higher, with pressures of the order of 100 atmos 
and higher, the emissivity approaches unity (¢€ > 0.95), while tor 
temperatures of the order of 9500°K and lower, with pressures of the 
order of 10 atmos and lower, the emissivity approaches zero (€<0.05). 
The variations of the emissivity between these approximate limits are 
shown graphically as functions of temperature and pressure with the 
geometrical depth set at 50cm. The variation of the emissivity with 


geometrical depth is also shown graphically at 12, 600°K and 20 atmos. 


ili 


TABLE OF CONTENTS 


Acknowledgements i 
Abstract ii 
Table of Contents dil 
Je Introduction i 
II. Evaluation of the Spectral Absorption Coefficient, Ps , 

and of Optical Density, X a: 

fe Calculation of Pu. 3 

E, Determination of Optical Density, X 8 
lil. Results 10 
References 12 
Tables 13 
Figures 15 
Appendix os 





i. INTRODUCTION 

‘& rocket utilizing nuclear energy with hydrogen as working 
fluid has been considered in the past. It is therefore of some interest 
to obtain the emissivity of hydrogen atoms as a function of temperature 
and pressure in order to assess the importance of radiant heat trans- 
fer for high-temperature systems utilizing hydrogen. If the tempera- 
tures are sufficiently high, the hydrogen will exist almost completely 
in the atomic state, and hydrogen atoms will make the rnajor 
contribution to the total emissivity. Furthermore, the treatment for 
hydrogen atoms provides a simple example of emissivity calculations 
for atomic systems. 

The emissivity, € , of a body is the ratio of the emitted radiant 
energy to the radiant energy emitted by a blackbody. The contributions 
to the emissivity of hydrogen atoms come from two sources: the con- 
tinuous spectrum and the discrete transitions. 


The blackbody curves may be obtained by plotting 


R° 2 27 Ae/us? 
w exp( Whew/kT) -1 against Ww for various temperatures. 


Here Ry dw (exb/tm" Bee) is the radiant energy emitted by a black- 
body, per unit area, per unit tire, in the wave number interval das (¢m') 
into a solid angle of 2MT steradians, T (°K) is the absolute temperature, 
C (cm/sec) is the velocity of light, By (erg-sec) is Planck's constant, 
and A (erg/°K) is the Boltzmann constant. If one considers the system 


of heated hydrogen atoms to be isothermal and uniform, then 


R=R., [t-expCR.X)] . This function can be plotted as a function of w on 


the same graphs as Ry. Here R,.dw (erg/em~-sec) is the radiant 








eee 


enerpy emitted by the hydrogen atoms, per unit area, per unit time, in 
the wave number interval dw, into a solid angle of 2 Tr steradians, 
-l er _— 
Pu (cm ~“-atmos ‘ is the spectral absorption coefficient, and 
4X (cm-atmos) is the optical density. The spectral absorption coefficient 
d »¢ d 
Fu is set equalto R,+R, , where Py arises from discrete 
, c ; saci 
transitions and PF, from continuous spectrum contributions. To 
obtain the emissivity one takes the ratio of the actual area under the 
Ruvs. W curve to the blackbody area, « sl » where 
-5 2 90.,4 
g- = 5.6699x10 ~ erg/cm”- K-sec. Examples of plots of Rw and 
Re, as functions of W are reproduced in Figs. 1 to 5 for various 
temperatures at a representative pressure of 40 atmos. For conven- 
° , | ° ° a 
ience of presentation, Ry/Rumax 22d Ru/Ru na, have been used as 
° © 
ordinates, where Ru,max is the maximum value of R,, . The con- 
siderations used for the calculation of FR, andhence R,, will now 


be described. 





Tl, EVALUATION OF THE SPFCTRYL ABSORPTION COEFFICIENT, 


Po , AND OF OPTICAL DENSITY, X. 


a. Ualculation ot P, 








"yy * bd * e * ra 
The mass absorption coeificient in cm” per gram of neutral 
hydrogen atoms for the continuous spectrum, X, , bas been given by 


Unsold! I) ag 


iL, 4b 








6+ m = c Le 
: cee ee iT a. + > (oe (1) 


lHiere wW 3g the mass of the electron (a), oy is the charge of the 
electron (esvu), mM, is the mass of the hydrogen atom (9); 
AA = new / RT » Up* Rhc/n'k is , Ris the Rydberg constant (cmon 
and y. is the principal quantum number of the hydrogen éetom. fo 
convert ae to P , which has the units of con” iteenee tie one must 
multiply by the atomic weight of hydrogen and divide by tne absolute 
temperature and the molar gas constant, R » in cin? el mg ; 
The contribution of the H'th spectral line whose center lies at 
We. can generally be determined by computing the rnaximuim value ot 
R,(K)=R,., (K) . which is usually the value of R. (K)at the line certer 
[= Ry. (4K) ] and multiplying this quantity by the intensity wicth defined 
as the wave number range for which R,(K) 2 4 Rwax(K) - Detailed 
calculations have shown that this procedure will give only an error of 
approximately one percent for the contribution of the line to the total 


emissivity. or collected methods of computing the spectral absorp- 


tion coefficient of isolated lines see, for example, “enner and 





ay F. 


Kavanagh!@). The maximum absorption coefficient for purely Doppler 


broadened lines, P , is given by the equation 


P = Coal Gaye’ / ane eye (2) 


-2 =) . ; 
where 5S (cm “-atmos 7 is the integrated intensity of the spectral 
line under consideration. The integrated intensity, S , may be 


written as 
S= SCN /N)[L- exe -Irewo/ RT) (3) 


where N is the total number of atoms per ee per atmos, Nv is the 
number of atoms in the lower state per oun” per atmos. Here S' is 


given by 


S=mTmaNf/mc* (4) 
where f is the "oscillator strength", which has been tabulated by 


Gaeui@ +), At equilibrium 


N./N = (4, /Q) exp (Ae ay) 1} (5) 


where a. is the degeneracy of the lower state involvec in the trans- 
ition, and Q is the partition function. In this cage 9, /Q is taken 

as n and D. at i/) where n is the principal quantum number 
of the lower state in the transition. In the present calculation the 
partition Ameen is, for practical purposes, equal to the contribution 
from the ground state. For the temperatures and pressures considered 
here it may be shown, using the estimates of Urey! >), that this 


procedure does not introduce errors greater than about one percent. It 





5. 
is now possible to write Eq. (3) as 


S= Sv? [1 a pi-hew./ fT )] exp ((1- 1f'X- R4Avc/RT)]. (6) 


To estimate the intensity width, 2|we, - We | , at half the value 
of R, 4. (K) , one uses the well known curves of growth (see, for 
example, Unsaid?) or Penner and Kavanagh! ey, Here w, isa 


wave number suchthat R,,(K\) = 4R_. (K) . The line shape para- 


max 


meter is given by 


a= 4.(In 2)‘/ brs, (7) 


where the natural half width has been neglected compared to the 


collision half width, &. The Doppler half width, Ry, , is given by 
a /2 
fy, = (2 kT Lv 2. /m,c*) Wo, (3) 


An estimate of the collision half width is obtained from the formula of 


Mivin' edo pt" ), 


= qe wft/2rme'we, (9) 


where W is the number of molecules per.cm?, Combining Eqs. (4) 


and (9), one obtains 
Ap.= TR /10 wo , ) (10) 


U 
where FP; is the pressure in atmos. Knowing a and PX, itis 


possible to read from the curves of growth the value of the quantity Le 


$= A(2n2)*/2 Ro b, (11) 





where 


A= 2] Wp - We| Rey. (12) 
Combining qs. (11) and (12), one obtains 


—fwr- we l= 5 2,/Unr Ne (13) 


If the Stark effect due to ionization is considered, P,,* may be 
taken as a first rough estimate, as Po + PR ,» where Pa is the 
absorption coetficient due to ionization broadening alone, and PR, is 
the absorption coefficient due to broadening without considering 
1onization. A relation for the absorption coefficient due to ionization 
broadening is given by 
oe KX (OA) 2” exp (hew. /RT) No (14) 

RT 
where N, is Avogardo's number, and, A(An”) (cm”/atom} is the 
absorption coefficient per absorbing atom at a wave length ar lA) 


from 1/W,. The relation for &% (4X ) is 
_ Ife! i 
A(AA) = 4.21 CuO N'/ (ary # (15) 


3 : 
Where WN , the number of protons per cm’, may be determined from 


. : 5 
the Saha equation, and Cy has been tabulated by Aller and 


* For amore ccornplete discussion of the Stark effect due to ionization 
see Ref. 5. Here broadening due to electrons is neglected. The 
electron effect has been discussed recently by B. Kivel, 5S. Bloom, 
and H. Margenau, Phys. Rev. 98, 495 (1955). 
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Unadia’?), Using Eqs. (14) and (i5), a check of ionization broadening 
has been made at various temperatures and pressures using Ha as 
anexample. In all cases it was foundthat [1-exe (Pox << 1 at 
the value of Wy previously determined without considering Stark 
broadening, so that it seems justifiable to take oe io except 
perhaps for the Lyman lines, where the factor exp (-hew, /RT) is 
unity. 

The method described above to determine contributions to the 
eniissivity due to discrete transitions does not apply in the case of the 
Lyman lines because they are too intense, and the blackbody curve 
along the intensity width of the lines varies appreciably in this area. 

At the temperatures considered here it is only necessary to determine 
modifications of the above method for Lyman 4& , since the succeed- 
ing Lyman lines give negligible contributions because of the small value 
of the blackbody radiancy, Ry , atlarge W . When ( W- Wz) is, 
for example, of such magnitude that it does not overlap an adjoining 
line, an expression for the discrete pressure broadening absorption 

(2) 


coefficient, P,; , has been given by Born For large values of 


&/a , Born's relation is 


m ay t-QMbY ge 1-(S)/ays (S) ele) 


4-(2)(6/a) + (48/2 (16) 
Be =(ae tas} ~ Hie) Mea +a lt + (§/ay')" ’ i 


4 

/2 
where § = [(w-w.)/ 4,1 (Ly 2) . When (W- W,) is, for example, of 
such magnitude that it overlaps an adjoining line, Born's method for 


pressure broadening and the treatment described above for the Stark 





ao. 


effect do not apply because one can no longer consider the contributions 
as resulting from "small" perturbations of well defined energy states. 
The problem of the very broad Lyman % line has caused difficulty in 
related astrophysical applications, and an exact treatment is not 
feasible at present. Therefore, a reasonable cut-off point has been 
chosen beyond which the above treatment is not justifiable. The cut-off 
point to the left of the line is taken arbitrarily as one half the distance 
between ihe Lyman 4 and Lyman § line centers. The cut-off point 
to the right of the line is unimportant since the blackbocy contribution 
is negligibly small. The application of Egs. (14), (15), and (16) to the 
Lyman lines results effectively in adding a small contribution of 
blackbody radiation from the left cut-off point to Wc , It is felt 
thai the neglect of the Lyman % contribution from ™=O to the left 
cut-off point does not give appreciable error, but this approximation 
may well deserve further investigation. 

The contributions due to the higher series, Brackett, Pfund, 
etc., have not been considered in the present computations. These 
series will give appreciable contributions only at much lower tempera- 
tures than were used. Even if these contributions were blackbody from 
the left of the first Paschen line to W=:=0 , the emissivity would in no 


case be increased by more than about Il °/0 ; 


B. Determination of Optical Density A 
The optical density, A , is defined as p= where p, is the 


partial pressure of the hydrogen atoms (atmos) and J! is the geometric 


depth (cm). The partial pressure is given as 





=O 


P= (24 /lt 2) Pe (17) 


where .X* is the degree of dissociation of hydrogen molecules and P, 
is the total pressure. The degree of dissociation, 2 , may be deter- 


mined from 


K.=@x/i-* )e, (18) 


where K, is the equilibrium constant for the dissociation of 
molecular hydrogen. Good estimates of K, can be obtained by ex- 
trapolating data given in the National Bureau of Standards Tebhep 
Values of * and f, are listed in Table I. From this table it is seen 
that hydrogen atoms may be expected to contribute significantly to the 
total emissivity. In a more exact treatment it would be necessary to 
consider the contributions of species other than hydrogen atoms which 


are present in relatively small amounts at the temperatures and 


pressures considered in the present investigations. 
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Wil, RESULTS 

By use of the above procedures the total emissivity of hydrogen 
atoms in equilibrium with hydrogen molecules at various temperatures 
and pressures has been computed. For most of the computations a 
characteristic geometrical depth of 50 cm has been used, and Figs. 1 
to 7 refer to that depth, as does the following discussion except where 
explicitly noted otherwise. In the ranges of ten.peratures and 
pressures chosen, the emissivity varies significantly with changes in 
the state variables. Typical plots used to determine the emissivity are 
shown in Figs. 1 to 5, where the outer curve represents the blackbody 
curve, the cross-hatched area represents contributions from the 
continuous spectrum, andthe solid area represents added contributions 
frorm the discrete transitions. The discontinuities in the continuous 
spectrum portions of Figs. 1 to 5 correspond to the onsets oi the 
various continua; the fund continuum begins to contribute at w = 4380, 
the Brackett continuum at W = 6850, the Paschen continuum at 
aw = 12190, andthe Balmer continuum at W = 27400. Itis not 
necessary to consider the Lyman continuum for the temperatures used 
here. As seen from the figures, the continuous spectrum contributions 
to the emissivity increase as the temperature increases. The signifi- 
cant lines are located as follows: Paschnen series, W,. = 5330, 7790, 
9130, 9940, and 10480; Balmer series, Wo= 15220, 20580, 23050, 
24400, 25150, and 25750; and Lyman series, W,= 682200. As may 
also be seen from the figures, the intensity width increases as the 


temperature increases, and hence discrete transition contributions also 





=P 


increase. In general, the continuous spectrum gives the major con- 
tribution to the emissivity at high temperatures (T > 920C°), and the 
discrete transitions give the major contribution at low temperatures 
(T=9200°K). In a similar manner it may be shown that for a constant 
temperature the emissivity increases with pressure because the 
number density of emitters increases (fixed 2, X increasing) linearly 
with pressure; furthermore, the lines are broadened as the pressure 
is raised. 

The data listed in Table Il are obtained from such plots as are 
shown in Figs. 1 to 5. These data are shown in graphical form in 
Fios. band 7. from these figures it is seen that the emissivity 
approaches unity (€ > 0.95) for temperatures of the order of 12, 500°K 
and higher, with pressures of the order of 100 atmos and higher, while 
for ternperatures oi the order of 9500°K and lower, with pressures of 
the order of 10 atmos and lower, the emissivity approaches zero 
(€ < 0.05). The graph in Fig. 8 shows the change in the total 
emissivity as a function of geometrical depth for constant temperature 
el = lize 600°) and pressure (P, = 20 atmos). If either the tempera- 
ture or pressure is decreased, the emissivity decreases and the 
é vs ! curve lies under the curve shown. Similarly, if either the 
temperature or pressure is increased, the «€ vs £ curve lies above 


the curve shown. 





(1) 


(2) 


(3) 
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TABLE II. EMISSIVITY OF ATOMIC HYDROGEN AS A FUNCTION 
OF TEMPERATURE AND TOTAL PRESSURE 


(L = 50 cm). 
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APPENDIX 


&. Sample Calculation 
A detailed computation is presented here in order to show an 
application of the previously described procedure more clearly. This 


specific calculation is for the emissivity at 10080°K and P, = 40 atmos 


with 2 = 50cm. The values used for the various constants are as 
‘follows: 
C= 2.9979 x eon cm/sec 
A= 6.6238 x co erp-sec 
A = 153803 Te we erg/°xK. 
ao = 5.6699 x 107° erg/cm“ssec-°K 
G = 4.8022 x 10 wag 
m= 9.1072 8 10°78" & 
: m,= 1,6819 x 10° = g 
R= 1.0968 x 10° em"! 


' 3 i 


R= 82.06 cm -atmos /mole-°X 


1. Blackbody Curves. The blackbody curves are plotted using 


the relation R23 =21h °° ie / fexp (hew/kT) - ic For T = 10080°K 


4 


this becomes Riu, = 3.74x 107° w °/exp (1.427x10° “w)- 1. Values of 


Ri, vs w are shown in Table A-I. The area under the blackbody 


1} 


curve, Se is 5.85x10 um cm” - bec. 


2. Contributions from the Continuous Spectrum. The rnass 





absorption coefficient [Eq. (1)] at 10080°K is given by 





aA 


x, = (462% 10*)[4- expCl427« 10%) G/w?® (sae 


=e Als (1 
where G-=ze /yve + @ 2.4, . Values of G obtained from Unsold 


are as follows:* 


w/10> G 
4.38 to 6.85 0.0596 
6.85 to 12.19 0.101 
12.19 to 27.4 0. 309 
27.4 to 109.7 6. 43 


Cc 
The spectral absorption coefficient, P, , is 


Po = X/, (-008)/R'T, (A-2) 


from which P, = 1. ZiGanor’ 


“y . The optical density, X = Pa LF 
is 1998 cm-atmos for this case ( fF, = 39.95 atmos). The term 

[1 - exp (- P, X)] can now be determined, and hence Ro can also 
be obtained. Table A-II lists values of Py and Ros for various W , 
The area of contributions from the continuous spectrum is found by 


graphical methods to be 1. 118x10/? erg/cm”-sec. 


3. Contributions from the Discrete Transitions. To calculate 


p' as given in Eq. (2), it is first necessary to compute 5 (Eq. (4)] 


* Unsold does not list values of G for 9200°K and 11300°K; however, 
these can easily be calculated as indicated by the formula for G. 
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and = S [f&q. (6)}. The value of N (iain w/ com g etme w) in Eq. (4) 
can be obtained from the perfect gas law, Prv= nRT, where 7n is 
the number of moles and Vv _ is the volume (cm?). Rewriting this 
with R= N. & » one obtains N= nN. /Fev = et . Hence Eq. (4) 


becomes 


S'= 645 «10° f, (A-3) 


where f is taken from the tabulated values of Unsaid"), The 
relations for S for the Balmer and Paschen lines, respectively, 
become from Eq. (6), 
S = 3.092% 10 “[1- exp(-6.427 * 15*w,)] S', 
and ; (A-4) 
S = 8.046% to °[4 - exp (-1.427 * 10 *w)] oe 


Equation (2) is then written 


==) Sloe 1o* S/We. (A-5) 


At the given temperature and pressure dry = 3. 58x107> We [Eq. (8)| ; 
PA 


° 


ana Q= 4.71x107 S / w (Eqs. (7), (8), and (9)} . Theterms Q 
and log 16 10.6 PX can now be determined and § _ is then found 


from the curves of growth. The intensity width from Eq. (13) is now 
2| Wy - Wo| = 2.403 § ks, (A-6) 


Table A-I]I gives the necessary data to obtain 2 |W, - aye | for the 
spectral lines considered, and the last column of this table gives the 


value of 2 | Cy, - We | computed from these data. The additional area 





es oe 


contributed by each line is usually obtained by multiplying 2| Wy - | 

by (Rs a R.,) at We. If 2| wp - W) o| is so large that either Ris 

or Rr deviate appreciably from linear variation along ne entire 
intensity width, then the area must be found by graphical methods. These 


additional areas are listed below: 


Balmer Paschen 
Line Area (erg/cm’-sec) Line Area (svg (em ee 

3 3. 1a gaio 4 0.635 x 10/9 

4 0.650 0.204 

5 0.224 0.125 

6 0.098 0.052 

i. 0.656 0.037 

8 0.037 


The Lyman contribution is taken as blackbody from the left cut-off 


point (W = 7. 46x10) to W=co, This area is found to be 


10 


0. 250x10 ame/ cm” 20 82: 


Hence, the total area contributing to the emissivity, i.e., the 


sum of the continuous spectrum and discrete transition contributions, 


10 1) 


is 16. 43x10 erg/cr“-sec, and therefore € = 1.643010" /S amma = 


0.281. 








227 = 


TABLE A-I, Ry AS A FUNCTION OF WwW ( T = 10080°K). 


w/10> Re/10° 
1 0.244 
3 1.90 
4,38 2.68 
6.85 7.30. 
10 11.80 
12.19 14, 40 
15 16. 83 
20 18, 22 
25 16.94 
27.4 15.70 
30 14,02 
32.5 12.44 
35 11.05 
40 7.87 
45 5.50 
50 B74 
60 1.55 
70 0.589 
80 0.212 
30 0.072 
100 0.025 


=—2o= 


TABLE A-ll. Ps AND Ra AS FUNCTION OF W (T = 10080°K, 


P,= 40 ATMOS, £ = 50 CM). 





W)/10° ps/107* Ro/10° 
4,38 3.82 1.94 
6.85 1.35 1.73 
6.85 2.29 2.68 
10 0.900 1.94 
12.19 0.540 1. 47 
12.19 1.65 4,05 
15 0.944 2.90 
20 0.426 1.49 
25 0.225 0.748 
27.4 0.172 0.531 
27.4 3.58 8.05 
30 2.75 5.94 
32.5 2.18 4, 40 
35 1.75 3,25 
40 1.17 1.65 
45 0.826 0.839 
50 0. 600 0.423 
60 0. 348 0.104 
70 0.220 0.025 
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